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Abstract 

We study a diffuse interface model for the flow of two viscous incompress- 
ible Newtonian fluids in a bounded domain. The fluids are assumed to be 
macroscopically immiscible, but a partial mixing in a small interfacial region 
is assumed in the model. Moreover, diffusion of both components is taken 
into account. In contrast to previous works, we study a model for the general 
case that the fluids have different densities due to Lowengrub and Truski- 
novski [27]. This leads to an inhomogeneous Navier-Stokes system coupled 
to a Cahn-Hilliard system, where the density of the mixture depends on the 
concentration, the velocity field is no longer divergence free, and the pressure 
enters the equation for the chemical potential. We prove existence of unique 
strong solutions for the non-stationary system for sufficiently small times. 
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1 Introduction 

In this article we consider a so-called diffuse interface model for two viscous, incom- 
pressible Newtonian fluids of different densities. In the model a partial mixing of 
the macroscopically immiscible fluids is considered and diffusion effects are taken 
into account. Such models have been successfully used to describe flows of two or 
more macroscopically fluids beyond the occurrence of topological singularities of the 
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separating interface (e.g. coalescence or formation of drops). We refer to Anderson 
and McFadden j8] for a review on that topic. 

The model which we are considering leads to the system 

pd t v + pv • Vv - div§(c,Dv) + Vp = - div(a(c) Vc <g> Vc) in Q T , (1.1) 

d tP + div(pv) = inQ T , (1.2) 

pd t c + pv • Vc = div(m(c)Vp) in Qt, (1-3) 

pfl = -p- 1 ^- (p + $( c ) + a(c) |7°H + 0(c) - a(c)i div(a(c)^Vc) (1.4) 



dc V v ' 2 

where = O x (0, T) and C M d , d = 2, 3, is a bounded domain with C 3 -boundary. 
Here v and p = p(c) are the (mean) velocity and the density of the mixture of the 
two fluids, p is the pressure, c is the difference of the mass concentrations of the two 
fluids, and p is the chemical potential associated to c. Moreover, 

S(c,Dv) = 2i/(c)Bv + 77(c) divvl, (1.5) 

where S(c, Dv) is the stress tensor, IEDv = |(Vv + Vv T ), v(c), 77(c) > are two 
viscosity coefficients, and m(c) > is a mobility coefficient. Furthermore, $(c) is the 
homogeneous free energy density for the mixture and 0(c) = $'( c )- 

This is a variant of the model proposed by Lowengrub and Truskinovski [27] for 
an interfacial energy of the form 
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f f |Vc| 2 

-Efree(c) = / $(c) dx + / «(c)— — (ix 



where the choice a(c) = p(c) was proposed in [27J. A derivation of the latter system 
can also be found in [U Chapter II] (in an even more general form). We note that 
the term a(c) Vc <g> Vc comes from an extra contribution to the stress tensor, which 
models capillary forces in an interfacial region. The model is a generalization of a 
well-known diffuse interface model in the case of matched densities which corresponds 
to the case p(c) = const., cf. e.g. Gurtin et al. [T9] . 
The system is equivalent to 

pd t v + pv ■ Vv — div S(c, Dv) + pV<7 = PPoVc in Q T , (1.6) 

c\p + div(pv) = mQ T , (1.7) 

pd t c + pv • Vc = div(m(c)Vpo) in Qt, (1.8) 

P^o + P 2 P = ^P 2 ^o-a(c)5div(a(c)5Vc) + 0(c) in Q T , (1.9) 

together with 

/ p (t) = / ^o(*) ^ = for all t e (0, T), (1.10) 
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where po is the mean- value free part of p, and p is a constant (depending on time), 
which is related to the mean values of the pressure and the chemical potential. Here 
p and g are related by 

$(c) a(c)\Vc\ 2 p 
9=— + L + - - /ic, 

and /i = po + p, p = J n p dx, g = go + g, g = f n g dx. Details on this equivalence 
can be found in [2J Section 3] . Here it is assumed that the fluids mix with zero excess 
volume, cf. [27]. This implies 

1 1+c + 1_c 



/3(c) 2pi 2p 2 
where pj is the specific density of fluid j = 1,2. Hence p(c) is of the form 

P(c) = — ^ (1.H) 
a + pc 

for a > and \/3\ < a. (More precisely, /3 = ^ - = ^ + 

We close the system by adding the boundary and initial conditions 

n-v| an = (n ■ S(c,Dv)) r + 7 (c)v T | sn = on S T , (1.12) 

d n c\ 9 n = d n p \ 9n = on^, (1.13) 

(v,c)|t= = (v ,co) in J), (1.14) 

where = dQ x (0,T). I.e., we assume that v satisfies Navier boundary condi- 
tions with some friction parameter 7: R — > [0, 00) and assume Neumann boundary 
conditions for c and p. 

In the case of matched densities, i.e., p = const., (3 = 0, resp., results on existence 
of weak solutions and well-posedness were obtained by Starovoitov [33J, Boyer [10J, 
Liu and Shen [26], and the author [4j. The long time dynamics was studied by Gal 
and Grasselli [HI HHJ HE], Zhao et al. [3E] and in [31 H]. Moreover, in [TT] Boyer 
considered a different diffuse interface model for fluids with non-matched densities. 
He proved existence of strong solutions, locally in time, and existence of global weak 
solutions if the densities of the fluids are sufficiently close. In the case of general 
densities, existence of weak solutions of a slightly modified system was shown in [2] , 
where the case of a free energy of the form 

_Ef ree (c) = / $(c) dx + / a(c) dx 

Jq Jq q 

with q > d is considered. This is the only analytic result for the model (ll.ip - (ll.4p so 
far known to the author. To the authors knowledge there are no numerical studies 
of this model in the case /3 7^ 0. A simplified model was used by Lee et al. [2] 
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in numerical simulations. Moreover, A. and Feireisl [6] constructed weak solutions 
globally in time for a corresponding diffuse interface model for compressible fluids. 

In the following, we will only consider the case p(c) ^ const., i.e., j3 ^ and will 
restrict ourselves for simplicity to the case a(c) = m(c) = 1. Moreover, we do the 
following assumption: 

Assumption 1.1 Let VL C M. d , d = 2,3, be a bounded domain with C l '-boundary, let 
e > 0, a > 0, and (3^0 such that \/3\ < a. Moreover, we assume v, 77,7 G C 2 (R) 
such that inf seM z/(s),inf seR ?7(s) > 0, inf seM 7(s) > 0, $ G C 3 (R) and p G C 3 (R) 
such that p(s) = if s G [— 1 — e , 1 + e ] for some e and p(s) > for all 

s G IL Finally, i/ inf seR 7(5) = 0, then we assume that Q has no axis of symmetry, 
cf. Appendix for details. 

Now our main result on short time existence of strong solutions is: 

THEOREM 1.2 Let v G i^(fi),c G H 2 (tt) with \c (x)\ < 1 almost everywhere 
and d n c \QQ = 0, d = 2,3, and let Assumption \l.l\ hold. Then there is some T > 
such that there is a unique solution v G if^O, T; L\ (O)) n L 2 (0,T; H 2 (Vl) d )), c G 
H 2 (0,T; H-^ity) C) L 2 (0,T; H 3 ({1)) solving (ir^ - (ID^ . p^ - (ITI^I . 

Precise definitions of the function spaces are given in Section [2] below. 

The theorem is proved by linearizing the system suitably, proving that the lin- 
earized operator defines an isomorphism between certain L 2 -Sobolev spaces, and 
applying a contraction mapping argument. To apply this general strategy it is essen- 
tial to reformulate the system (ll.6l) -( fL9l) . To this end, we eliminate po an d go first 
from the system. Then the principal part of the linearized system (around (vo,co)) 
is 

a t v-div§(c ,Dv) + -^-Vdiv(p 4 Vc') = fi inQ T , (1.15) 

pa 

^c'-^divv = h inQT, (1.16) 

where d ~ pc and p — p(co), cf. Section |3] below. One of the essential steps in the 
proof of the main result is the analysis of this linearized system. To this end we split 
v in a divergence free part w = P CT v, a gradient part VG(div v), which is determined 
uniquely by g = divv, and a lower order part, cf. Section H] for details. A crucial 
observation is that c', which is related to divv via fl!.16j> . solves a kind of damped 
wave/plate equation. More precisely, c' solves an equation of the form 

d 2 c' - A(a(c Q )d t c') + div(po 4 Vc') = / (1.17) 

up to to lower order terms for some a(co) > 0. In order to solve this equation we 
will apply the abstract result of Chen and Triggiani [12]. - We note that the same 
kind of linearized system arises in Kotschote [22] , where existence of strong solutions 
locally in time is proved for a compressible Navier-Stokes-Korteweg system. 



5 



Remark 1.3 It is interesting to compare fll.15p - fll.16p to the linearized system of 
the Model H for the case of matched densities, i.e, (ll.ip - (ll.4p in the case when (3 = 
and thus p(c) = const. Then the pressure p is no longer part of f ll.4p . divv = 0, p 
can no longer be eliminated from the system, and the principal part of the linearized 
system is 

d t Y - div§(c , Dv) + Vp = fi in Q T , 
divv = in Q T) 
d t c + div(m(c ) VAc) = / 2 in Q T . 

Hence the linearized system is very different. In particular, the principal part for c is 
given by a fourth order diffusion equation with a(— div(m(c)VAc)) C (—00, 0] for a 
suitable realization. While the corresponding operator A\ to (11.171) (after reduction 
to a first order system, cf. Section |4|) still generates an analytic semigroup, but the 
spectral angle 5 < | can be arbitrarily close to | in certain situations, cf. Remark 14.31 
below. Moreover, note that the Cahn-Hilliard part is decoupled from the Navier- 
Stokes part on the level of the principal part of the linearized system in the case 
= 0, which is no longer the case if (3 7^ 0. - We hope that the insight on the analytic 
structure of the system fll.6p - fll.9l) will help to create stable numerical algorithms, 
which are not available so far to the best of the author's knowledge. 

The structure of the article is as follows: In Section [2] we summarize some notation 
and preliminary results. The main part of the article consists of Sections [3] and HI 
First, in Section [3] the system is linearized and the contraction mapping principle 
is applied on the basis of the well-posedness result Theorem 13.11 for the linearized 
system. Afterwards, in Section @] this result is proved. 

2 Preliminaries 

Notation: Let us fix some notation first. For a, b G W 1 let a g) b G W. dxd be defined 
by (a ® b) itj = aibj. Moreover, for A, B 6 R dxd let A : B = ti(A T B) = £ij=i a^hj. 
In the following n will denote the exterior normal at the boundary of a sufficiently 
smooth domain Q C lR d . Furthermore, f„ := n ■ f and f r := (J — n (g) n)f = 
f — f„n denote the normal and tangential component of a vector field f : dQ — > M. d , 
respectively. Furthermore d n := n-V, V T := (/ — n(g>n)V, and d T . := ey (I— n(g)n)V, 
j = 1, . . . , d, where ej denotes the j-th canonical unit vector in M. d . If v G C 1 (r2) d , 
then Vv = Dv = (d Xj Vi)fj =1 denotes its Jacobian. Moreover, if A = (ay)^ =1 : Q — > 
R dxd is differentiable, then div A(x) := ($^j=i d Xj aij(x)) d =1 for all x G Vt. 
If X is a Banach space and X' is its dual, then 

(f,g) = (f,g)x>,x = f(g), fex', g ex, 

denotes the duality product. The inner product on a Hilbert space H is denoted by 
(., .)h- Moreover, we use the abbreviation (., .)m = (•, -)l 2 (m)- 
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Function spaces: If M C R d is measurable, L q (M), 1 < g < oo denotes the usual 
Lebesgue-space and ||.|| 9 its norm. Moreover, L q (M;X) denotes its vector-valued 
variant of strongly measurable g-integrable functions/essentially bounded functions, 
where X is a Banach space. If M = (a,b), we write for simplicity L q (a,b;X) and 
L q (a,b). 

Let Q C R d be a domain. Then W™(Q), m G No, 1 < g < oo, denotes the usual 
L^-Sobolev space, W™ (Q) the closure of C °°(Q) in W 9 m (ft), W" m (fi) = (W£ (fi))', 
and W"- m (n) = (Wpln))'. The L 2 -Bessel potential spaces are denoted by H s (tt), 
sGK, which are defined by restriction of distributions in H s (R d ) to Q, cf. Triebel [3"51 
Section 4.2.1]. Moreover, H S (M) denotes the corresponding space on a sufficiently 
smooth compact manifold M. We note that, if Q C R d is a bounded domain with 
C 0,1 -boundary, then there is an extension operator Eq which is a bounded linear 
operator E n : W™{Vl) -> jy p m (R d ), 1 < p < oo for all m G N and E n f\ n = f for all 
/ G W^ 71 ^), cf. Stein [MJ Chapter VI, Section 3.2]. This extension operator extends 
to E n : H s (tt) -»■ # s (R n ), which shows that # s (fi) is a retract of # s (fi). Therefore 
all results on interpolation spaces of H s (R n ) carry over to H S (Q). We refer to Bergh 
and Lofstrom [9] for basic results in interpolation theory. In the following (., .)[#] 
and (■, .)e, q w iU denote the complex and real interpolation function, respectively. In 
particular, we note that 

= (H S0 {tt),H s i(ty) e!2 = H s (n) (2.1) 

for all 6 G (0, 1) where s = (1 - 9)s + 0si, s , s x G R. 
Moreover, we define 

^(fi) = {uGi/ 1 (fi) (i :n-u|^ = 0}. 

The usual Besov spaces on a domain or a sufficiently smooth manifold are denoted 
by Bp q (Q), Bp (M), resp., where s G R, 1 < p, g < oo. For the convenience of the 
reader we recall that Bf^M = # s (R d ) for all s G R and 

B s p %(R d ) ^ B s Mi (R d ) ^ B; m (R d ) for all 1 < q x < q 2 < oo. 

Moreover, we have the Sobolev type embeddings 

fl£, ff (R n ) ^ ^(R") if si > ^0 and fll - i > a „ - 
B%?(R n ) ^ C b °(R n ) 

for any 1 < p, q < oo. Finally, due to Hanouzet [20j Theoreme 3] we have the useful 
product estimate 

\\fg\\m(R d ) < Cp||/|| s d/ P(K(i) ||g||_H-i( K d) (2.2) 

d 

for all / G -Bp ;1 (R d ), # G H 1 ^*) provided that 2 < p < oo, see also (2D Theorem 6.6]. 
All these results carry over to sufficiently smooth domains and d- dimensional mani- 
folds. 
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Let / = [0,T] with < T < oo or let I = [0, oo) and let X be a Banach 
space. Then BUC(I; X) is the Banach space of all bounded and uniformly continuous 
/: I — > X equipped with the supremum norm. The space of all (uniformly) Holder 
continuous functions /:/—>■ X of degree s G (0,1) is denoted by C s ([0,T];X) 
normed in the standard way. Furthermore, we have the useful embedding 

BUC([0, T};X 1 ) n C s ([0, T]\ X ) ^ C^([0, T]; X), 

where < s,0 < 1 provided that \\f\\ x < C||/||^||/||^ for all / G Xq fl Xl 

Finally, / G W/ fc (0,T;X), 1 < p < oo, A; G N„, if and only if g£ G 

L p (0, T; X), where -j* denotes the fo-th X-valued distributional derivative of /. Fur- 
thermore, we set H 1 ^, T; X) = W%(0, T; X) and for s G (0, 1) we define H a (0, T; X) = 
S| 2 (0, T; X), where / G B| >2 (0, T; X) if and only if / G L 2 (0, T; X) and 

11/11^(0,7^ = ii/n W) + [ jf u ®s$F X dtdT < °°- 

In the following we will use that 

io /o |t-r| 2 ^ 

< T f \t-r\ 2{s '- s) ^dtdr\\f\\ csl{W) < C s ,, s T^'^ +1 \\f\\ csl{W) 
Jo Jo 

for all < s < s' < 1, which implies 

II/||h.(owo < ^'^II/IIgV([o,t] ; x) for all / G C s '([0,T];X) (2.3) 
provided that < s < s' < 1. Finally, we set for s G (0, 1) 

Hi' s {S T ) = L 2 (0,T; H s (dQ)) fl Hi(0, T; £ 2 (<9ft)), 

where St = <9fi x (0,T) and Q is a bounded domain with C 1 -boundary. 
Now let Xo,Xi be Banach spaces such that X\ <-» Xo densely. Then 

Wp(I; X ) n U\l- X) ^ SC/C(J; (X , X)^^), 1 < p < oo, (2.4) 

continuously for / = [0,T], < T < oo, and / = [0, oo), cf. Amann [7J Chapter III, 
Theorem 4.10.2]. 

In order to solve the linearized system in the following, we will use the following 
abstract result: 

THEOREM 2.1 Let A: T>(A) C H — > H be a generator of a bounded analytic semi- 
group on a Hilbert space H and let 1 < q < oo. Then for every f G L q (0, oo; H) and 
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u G (H, V(A)) 1 _i there is a unique u: [0, oo) — > H such that ^jt, Au G L g (0, oo; H) 
solving 



du 
~dt 



(t) + Au(t) = f(t) for all t > 0, 



u(0) = u . 

Moreover, there is a constant C q > independent of f and u such that 
du 



dt 



+ \\Au \\L1 



Li{0,oo;H) 



(0,oo;H) 



<C n 



Li(0,oo;H) + \\ u 0\\ (H,V(A)) 1 _ 



Proof: In the case Uq = the statement is the main result of [13]. The general case 
can be easily reduced to the case uq = by subtracting a suitable extension. The 
existence of such an extension follows e.g. from [3 Chapter III, Theorem 4.10.2]. ■ 



Weak Neumann Laplace equation: In the following we assume that Q C R d is 
a bounded domain with C ' ^boundary. Given / G L 1 (r2), we denote by m(f) = 
W\ In f( x ^ d% hs mean value. Moreover, for m 6 1 we set 

L« m) (fi) := {/ G L«(fi) : m(f) = m}, 1 < q < oo. 

Then 

Pof :=/-m(/) = /-p^/(x) dx 
is the orthogonal projection onto LjL(f2). Furthermore, we define 

H{ 0) = Hfa(fl) = H\Q) n Lj^fi), (c,d) Hlo){n) := (Vc, Vd) L2(n) . 
Then ffL(fi) is a Hilbert space due to Poincare's inequality 

\\f-m(f)\\ LPi n)<C p \\Vf\\ LP{ n), 

where 1 < p < oo. Moreover, let HZs = HZ}(Q) = H^(£l)'. Then the weak 
Neumann- Laplace operator A^: HhJQ) — ► JT^ (fl) is defined by 

- (A N u, ip) h i q) H -i = ( Vw, Vip) for all if G (fi) . 

By the Lemma of Lax-Milgram, for every / G Hz}(ft) there is a unique u G H^(fl) 
such that — Ajyu = /. More precisely, — Ajv coincides with the Riesz isomorphism 
K: H^V) ^ HZ,](n) given by 

(ftc, d) H -i H i = (c, dW = (Vc, Vd) L 2, c, d G Hf 0) (n). 
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We equip H, Q JQ) with the inner product 

(f,g) H -i = (VA-V,VA-^) i2 = (A-V,A-^) H(lo) . (2.5) 

In particular this implies the useful identity 

((-A N )f,g) H -i = (f,g) L2 for all / G (fl) , g G L 2 (0) (fl). (2.6) 

Moreover, we embed (f2) and L 2 j (fi) into if™ (fl) in the standard way by defining 

( c > V) H 7 n \,H} m = / c(x)tp(x) dx for all y> G 
for c G fJL(O). This implies the useful interpolation inequality 

\\f\\ 2 L2 = -(VA N l f, V/) L2 < ||/||^||/||^ o) for all / G H^fl). (2.7) 
Furthermore, if A^u = f for some / G if™ (f2), then 

Nlnfodi) < \\f\\n- ]{ ny (2-8) 

We note that, if u G H^{fl) solves A N u = f for some / G Ll-j(fi), 1 < q < oo, 
and is of class C 2 , then it follows from standard elliptic theory that u G W 2 (fl), 
Au = f a.e. in fl, and <9 n ti|,9Q = in the sense of traces. If additionally / G Wg(fl) 
and dfl G C 3 , then u G W 3 (f2). Moreover, 

IMU^ ( n) < C q \\f\\ wm for all / G W/ fc (fi) n = 0, 1, (2.9) 

with a constant C q depending only on 1 < q < oo, d, k, and fl. 
Finally, we define div„: L 2 (fl) -» Hz}(fl) by 

(div n /, ¥>>H-i(n),Hi 0) (n) = ~(/> V^)L2 (n) for all G 
Note that A^w = div n Vw for all u G H^(fl). 

Helmholtz decomposition: Recall that we have the orthogonal decomposition 

L 2 (fl) d = Ll{Q)®G 2 {Q) 
G 2 (fl) = {V P eL 2 (fl):peHf 0) (fl)}. 

Here L 2 a (fl) is the closure of {u G C^(fl) d : divu = 0} in L 2 {fl) d . The Helmholtz 
projection P a is the orthogonal projection onto L 2 (fl). We note that P a f = f — Vp, 
where p G HhJfl) is the solution of the weak Neumann problem 

(Vp, V<p)n = (/, Vy>) for all ^ G C$(fi). (2.10) 

We refer to Simader and Sohr [31] and Sohr [321 Chapter II, Section 2.5] for details. 

We conclude this section with two technical results related to the Navier boundary 
condition (I1.12p . which will be needed in Section HI 
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Lemma 2.2 Let Q C M. d be a bounded domain with C k -boundary, k > 2. Then there 
is a first order tangential differential operator A = Ylj=i o J j{x)d Tj , aj G C fc_1 (9fi), 
such that 

(n • V 2 u) T \ dn = V T 7i« + A-f u for allu e H 2 (VL), (2.11) 
where 'jjU = d 3 n u\do,. 

Proof: Since V r = (I — n (g) n)V, we obtain 

d 

{n-V 2 u) T \ dn = (I-n®n)(n-V 2 M)|^ = (/-n®n)V(a n M)-^(a 7 - j n)-VM 

3=1 

for all u G H 2 (Q). Since d Tj n ■ n = |(9 Tj |n| 2 = 0, d T .n G C* -1 (<9JT) is tangential. 
Therefore ( 12. lip is valid. ■ 



Lemma 2.3 Let Q C M. d be a bounded domain with C 2 -boundary, < T < oo, 
v G C\R) with inf seM v{s) > 0, and c G H 2 {Vt), d = 2, 3. 

1. There is a bounded linear operator E: H2(dn) d -> # 2 (fi) d swc/j t/iat 

(n-2©£a) r | 9Q = a r , £a| sn = 0, div£a = 

for alia G H^(dQ) d . Moreover, there is a constantC > suc/i t/iat U-EaHij-im) < 

There is a bounded linear operator 

e t : m^(s T ) d -»■ ^(o.ri^^fniJ^o.riL^o)) 

(n • 2/y(co)©^ T a) r | an = a T , S T a| 9n = 0, div£ T a = 0, £ T a| i=0 = 

/or all a G H^'^(ST) d - Moreover, the operator norm of Et can be estimated 
independently of < T < oo. 

Proof: To prove the first part let A = Ea G H 2 {yt) d such that A\qq = 0, d n A\ d Q = 
a T and ||A|| H 2 (n) < C\\a\\ H x {gn) , \\A\\ H i m < C\\a\\ H _i {m) for all a G #§(<9ft) d . If 

Q = W 1 ^ 1 x (0, oo), the existence of such an extension operator E follows e.g. from 
McLean [281 Lemma 3.36]. From this the result for a general bounded C 2 -domain 
follows by standard localization techniques. 
Then we have 

(n ■ 23A) T \ dn = (V T A n + d n A T )\ dn = + a T , 
divA| 9Q = (div T A + d n A n ) \qq = 0. 
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Since div A|an = 0, div A G Hq (fl) PlL^ (Q) and we can apply the Bogovski-Operator 

B, cf. e.g. [17], to div A. Hence we obtain ^(div A) G H$(Q), divi?(div A) = div A, 
and 

^(divA)!!^) < C||A||^ ( n) < c 'M H h m y 
\\B(dw A)\\ HHn) < C7||A||fli ( n) < C'M H _ km) 

for all a G H^(dQ). Therefore A := A — B(div A) has the properties stated above. 

Finally, because of [5J Lemma 2.4], for every a G H*>2(S T ) d there is some A G 
£ 2 (0,T;# 1 (fi)) d nif 1 (0,T;£ 2 (fi)) < ' such that (n • i/(c )DA) T | an = a T , A| t=0 = 0, 
div A|^q = 0, Al^ = 0. Moreover, the extension can be chosen such that 

||A||i2( 0) T;H 2 ) + IIAH^i^l 2 ) < C , ||a||^i i i (s ,^ ) 

with C independent of T and a. Analogously to the first part div A G L 2 (0, T; HQ(Q)) d . 
Hence E T a : = A - £(div A) G L 2 (0, T; H*(tt)) d . Moreover, due to [HI Theorem 2.5] 
we also have 

||5(divA)||#x( 0) T;i,2(n)) < C\\ div n A|| Hl(0)T .„-i (Q)) < C'\\A\\ H i (0!T . L 2 in)) , 
where C, C > are independent of T. Altogether £Jy has the stated properties. ■ 



3 Short Time Existence of Strong Solutions 

In this section we prove existence of strong and unique solutions of system fll.6l) - fll.lOI) 
and (Il.i2p - fll.i4p locally in time in the case a(c) = m(c) = 1, i.e., prove Theorem 1 1.2 1 
As noted before we will assume that f3 7^ since in the case (3 = the linearized 
system is completely different and short time existence of strong solutions is known 
in that case, cf. e.g. [I]. In this case we can eliminate the generalized pressure go 
and the chemical potential p as follows: 

First of all, because of ( 11. lip , one easily calculates that 

d P o 2 9{pc) dp 2 

^ = -^' -d^ = P+ d-c C = ap - (3 ' 1} 

For the following let c be a sufficiently smooth solution such that \c(t, x)\ < 1 + Eo, 
where Eq > is as in Assumption 11.11 Then ( 11 .Th and ( 13.11) imply 

-(3p 2 (d t c + v- Vc) = -pdivv. (3.2) 

Combining this with ( II. 8ft . we obtain the simple identity 



div v = /?A/i . 



(3.3) 
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3 SHORT TIME EXISTENCE OF STRONG SOLUTIONS 



Thus fio = (3 1 A 7V 1 divv = {3 1 G(divv) since n ■ Vpolan = 0, where G(g) is defined 
by 



AG(g) = g in Q, 
d n G(g) = on ffl, 

and j n G(g) dx = 0. Note that this implies 

VG(divv) = (7 — P a )v. 



(3.4) 
(3.5) 



(3.6) 



Using this and (11.91) . we can eliminate go, /io and pit) from the system fll.6p - fll.9l) 
and obtain the equivalent system 

d t v + v • Vv - p- 1 div S(c, Dv) 

+ ^V(p- 2 (Ac-0(c))) =^G(divv)Vc-^V(p- 1 G'(divv)) in Q T , (3.7) 



ft 



together with 



pd t c + pv • Vc = f3~ 1 div v. 



ft 2 



mQ T , (3.8) 



n ■ v\ 9n = (n ■ S(c, Dv)) T + 7(c)v 



-nan 



9„.c| 



an 



(v ,c J 



on 5^, 
in £7. 



(3.9) 
(3.10) 



This is indeed an equivalent system since, if (v, c) solve the system above, we can 
simply define go and pit) by the equation (11.91) and /io by fj, = /3 _1 G(divv). Then 
n • Vpolan is automatically satisfied. 

We will construct strong solutions by linearizing the system, proving that the 
associated linear operator is an isomorphism between suitable L 2 -Sobolev spaces, 
and applying the contraction mapping principle to prove existence and uniqueness 
of the full system for sufficiently small times. 

To this end, let c G H\0,T ; H\Q)) n L 2 (0, T ; H 3 {Q) n H 2 N {Q)) be such that 
Co\t=o — Co- The existence of such an cq follows from Theorem 13.11 below. Then 
(I3.7p - ()3.10p are equivalent to 



Lie) 



7"(v,c), 



K p(c-c ) / 

where for given c the linear operator L{c): Xj- — >■ is defined by 



(3.11) 



He) { 1 



(d t v - div S(c, Dv) + ± V div(p- 4 Vc')\ 
9jc' — div v 
(n-S(c,D(P (T v))) r + 7( C )(P (J v). 

V (v,c)U / 



e x 7 



13 



and d corresponds to pc. Here we have used V(pc) = ap 2 Vc and set §(c, Dv) = 
2£(c)Dv + 77(c) divvl, £(c) = p(c) _1 z/(c), 77(c) = p(c) _1 ?7(c), and 7(c) = p(c) _1 7(c). 
Moreover, J 7 : Xt — > Yt is a non-linear mapping defined by 



^(v,c) 



*i(v,c) 



/ ^i(v,c) \ 

-pv • Vc - p<9 t c + <9 t p(c - c ) 

n • S(c, V 2 G(div v))) T + 7(c) V T G(div v ' 

OS. 

V (v ,o) y 

(c) G(divv) . _ 2 



9fl 



G(divv)Vc 1 
— h -V 

P 



P 2 P 



[p ,div]Vc 



+ -V(p~ 2 Ac ) - v • Vv - Vp- 1 ■ S(c, Dv), 
P 

and X T = X T x X T , 

X\ = {uE H 1 (0,T;L 2 ({l) d ) nL 2 (0,T;H 2 (ty d ):n-u\ dn = 0} 
{c G ^(0, T; ff 1 ^)) n L 2 (0, T; # 3 (fi)) 



c t=o 



0, n • Vc' 



an 



0} 



F T = L 2 (Q T ) d xL 2 (0,T;i7 1 (fi))x{aG J ?/i5( > S T ):a n = 0}x^(r2)x^(r2), 

where H l ^(S T ) := H^(0,T; L 2 (dty) n L 2 (0, T; H^(d£l)) Here [A, B] = AB - 
denotes the commutator of two operators. The spaces X T , X T , and Yt are normed 
by 



\ c \\x* 



L 2 (Q T ) 



||(f,#,a, V )||y T 



(^v,V 2 v)| 
\(c',d t c',d t Vc',V 3 c') 
( f , Vg)\\ L t { Q T) + ||a|| 



|| v |t=o||j3' 1 (n)) 

| l2(Qt) + ||^|t=o||ip(n), 



Hhh(S T ) 



|Vo||i?i(n) + || c 1| ^2(0) ■ 



In order to apply the contraction mapping principle to (13.111) for sufficiently small 
T > 0, it is essential that L(cq) is an isomorphism: 

THEOREM 3.1 Let c G H 2 (Q), let T > 0, and let Assumption O hold true. 
Then L(cq) : Xt —> Yt is an isomorphism for every < T < Tq and there is a 
constant C(Tq) > such that 

-ii 



\\L(c ) 



\C(Y T ,X T ) 



< C{T ) for allO<T< T . 



(3.12) 



The proof of this theorem is postponed to Section HI The second ingredient for the 
application of the contraction mapping to (13. lip is the fact that J 7 : Xt — > Yt is 
locally Lipschitz continuous with arbitrarily small Lipschitz constant if T > is 
sufficiently small: 

Proposition 3.2 Let R > and let Assumption be satisfied. Then there is a 
constant C(T,R) > such that 



||7-(vi, ci) - -F(v 2 , c 2 ) \\ Yt < C(T, R) || (vi - v a , d - c 2 ) \\ Xt 
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3 SHORT TIME EXISTENCE OF STRONG SOLUTIONS 



for all (vj,Cj) G Xt with || (v.,-, Cj) \\x T < R and Cj\ t =o = cq, where j = 1,2. Moreover, 
C(T,R) -> as T -> 0. 

Proof: Let F 2 (v,c) = — pv ■ Vc and let F 3 (v,c) = —d t p(c — c ). For F'(v,c) = 
(Fi(v, c), F 2 (v, c)) we will show that 

||^'(Vl,Ci)-F'(v2,C2)|| LP ( 0) T;i2(n)<«xJfl(«)) - C'CPj ^ T o) II ( v i ~ v 2 , c i - c 2) \\x T (3.13) 

for all < T < To and for some p > 2. (Note that the third component of J 7 is 
constant.) Then the statement of the proposition for these terms follows from the 
estimate 



i i 



L 2 (0,T;X) <T 2 v\\J \\LP(0,T;X)j 

where X is an arbitrary Banach space. In order to estimate the terms involving c, 
we use that 

\\c\\L°°(0,T;H 2 (n)) < C\\c\\ X 2 (3.14) 

due to (12 .4p and (12. ip . where C is independent of T > 0. Since H 2 (fl) is an algebra 
with respect to pointwise multiplication, we have F(c) G L°°(0, T; H 2 (Q))) for all 
F G C 3 (M), c G Xf,, as well as 

||F( Cl ) -F(c 2 ) \\ L ~ (0 ,T;H> m < C{R,F)\\c x -c 2 \\ xl (3.15) 

for all Cj G X\ with ||cj|| X 2 j = 1,2. Hence 

||V(/3(ci)" 2 0(ci) -/3(c 2 )~ 2 (/)(c2))|| L ^ ( o,T;Hi(n)) < C(R, (j), p) \\c x - c 2 \\x* 
for all Cj G X\ with ||cj|| X 2 < R, j = 1,2. Moreover, 

||(7(div v) 1 1 l°°(o,T;H2 (n)) < C|| vIU^o.tw 1 ^)) < C"ll v llxi 

for constants C, C" independent of T > due to (12. 9p and (12.41) . Since the product of 
Lipschitz continuous functions is again Lipschitz continuous, it is sufficient to verify 
that all the products appearing in F(v, c) are well-defined, which is done as follows: 

||G(divv)Vc|| L oo ((Vr . L 2 (f7)) < C||G(divv)|| L oo (0)T . H 2 (n)) ||c|| L oo (or . H 2 (n)) < C(R) 
||V(/o" 1 G ! (divv))||£<»(o 1 r;La(n)) < C||G(divv)|| I( <» ( o,T;tf«(n))||p -1 |U«»(o ) r;fla(n)) <C{R) 

||V • Vv||x,4(o,T;L2(fi)) < || v ||L°°(0,T;L6(n))||Vv|| L 4( 0iT;i 3(n)) 

1 1 

< C'(^)||v||^ oo(0)T ._ ff i (n)) ||Vv||| f2(0iT . I(6(n)) < C (R) 

II Vp _1 ■ S(c,Dv)|| L 4 (0 ,T;L2(n)) < C'||Vp" 1 || L oo ( o i T;L6)||(//(c),77(c))|| L oo ( Q T) ||Vv|| L 4 ( o,T;L3) 

< C^llVvll^^^llV^lli^^ < C\R) 

||pv • Vc||i4(Q )T;fl i(n)) < C|| v llL 4 (o,T;B.i a (n))ll Vr(c)|| L oo (0iT . H i) 

i i 
^ ^ll v llla( ,r;fla)||v||2 00(()ir . H ij||Vr(c)||i ( oo(o I Ti£ri) 

< C(R), 
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for all (v, c) G Xt with ||(v, c)||x T < -R, where r'(c) = p(c) and we have used ( 12 ,2p 
and 3 

In order to estimate V([p~ 2 , div]Vc), we use that V[p~ 2 , div] Vc = — V(Vp -2 ■ Vc) = 
V 2 p -2 • Vc + Vp~ 2 ■ V 2 c and 

l|VV 2 -Vc|| L4(0 

< C(i2) || Vc|| La( o, Ti Ha(n)) || Vc|| L co (0 ,r ; Hi(n)) < C'(R) 
i i 

due to ||/||z,°°(n) < C||/|| 3 < C'H/Hlfill/11^2- The same estimate holds true for 

B 21 (Q) 

Vp~ 2 ■ V 2 c. Hence ( 13131) holds with p = 4. 
Furthermore, we estimate i^v, c) as follows 

ll(p(Cl) - p(c2))^Co|| L 2 (0iT;H i {n)) 

< C||p( Cl ) -p(c 2 ) ||, 

II W(ci) - 9 t p(c 2 ))(c! - c )\\ L 2(o,T;m(n)) 

< C||a t p(ci) - dtp(c 2 ) || mo,T-m) || ci - ^11 £00(0,^/2) < C(i2,oo)T»||ci - c 2 || X 2 

||<9 4 p(c 2 )(Ci - C 2 )|| L 2 (0iT ;Hi(Q)) 

< C||^p(ca)||i2( 0j T;tfi(n))||ci - calli-^r;^/ 3 ^)) - c '(- R J Co) rs ll c i - C 2|U|, 
since 

x 2 C3([o,r];£r 1 (n))nL oo (o,T; J fr 2 (n)) c*([o,t]; (3.16) 

ifi(fi) M> and (ci - c 2 )| t=0 = (ci - c )| t=0 = 0. 

Finally, it remains to estimate the third component of J-(y,c). To this end we 
use that 

||(n.V 2 G(divv) T || L2(0iT;H 3 ( ^ )) < C||v|| L 2 (0iT;H 2 (n)) 
||(n-V 2 G(divv) T || H1(0iT;H _i (9Q)) < C||v||ff 1( o,T;^(Q)) 

for all v G X\ since (n • V 2 G(divv)) T = Ajqv for some first order tangential differ- 
ential operator A, cf. Lemma 12.21 Hence 

||(n.V 2 G(divv)) r || L2(o ^, (m)) 

< ^||(n.V 2 G(divv)|| BW];H , (an)) <CTh,\\ x , 

for all v G X^, due to (j2~3|) and 

|| (n- V 2 G(divv) T || i „ , < T*||(n- V 2 G(divv)|| i, r , „ „ < CT3||v|| x i 
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3 SHORT TIME EXISTENCE OF STRONG SOLUTIONS 



due to O and since BUC([0, T); H 1 ) fl C^([0, T); H~^) ^ Ci([0,T];L 2 ) because 



of \\f\\v < CH/ll^ll/llj * Together we have 

||(n.V 2 G(divv) r || //i , (5T) <C^||v|| x , 
By the same arguments one shows that 

||(VG(divv)) T || Hi , 4(fti) <CTi||v||^. 

Moreover, using (I3J6J), \\fg\\ H l {m) < ^ll/IU^^IM!^^ as well as \\fg\\ L \an) < 
||/IU*(an)||»IU*(«i) < C|l/ll^ (an) ll»ll^(an)' one obtains 
\mci)-Hc2))a\\ Hk{0T , L2{m)) 
< CT^\\{v{ Cl ) - ^(c 2 ))|| c i ([0jT] .^ (an)) ||a|| i2(0T;ff i ( ^ )) 
+CT 1 q(v(c 1 ) - ^(c2))|| ci([0iT];Bli(an)) ||a|| H|(0!T;i . 2(an)) 

and therefore 

||(P(ci) - % 2 ))a| 



, i 



< C^^T^Id -C2|| sac([0ir];H2(n))nC i ([0)T];Hl(n)) ||a|| H i,i (ST) (3.17) 

for all a G H^'^(St) and Cj G Xf> with ||c 3 -||x 2 < i? and Cj | i= o = Co, j = 1,2. 
Combining these estimates, we obtain 

||(/>(c 1 )n-V 2 G(divv 1 )) r -(z>(c 2 )n.V 2 G(divv 2 )) T || H x,i (5T) 

< ||((P( Cl ) - z>(c 2 ))n ■ V 2 G(div Vl )) r ||^,i (ST) 

+\\(u(c 2 )n • V 2 G(div( Vl - v 2 )) t || hH(St) 



< C(i2)T* ( || vi-v 2 1| x x + || ci-ca | 



for all Uj = (vj,Cj) G with ||Mj||x T < -R and Cj| t= o = c , j = 1,2 since 
||(Vv, v)||^i 1^ ^ < C||v||jjfi. In the same way one can estimate 7(ci)(VG(divv!) r — 

7(c 2 )(VG(div v 2 ) T |gn, which proves the necessary estimate of the third component 
of T. 

Altogether this proves the proposition. ■ 

Combining Theorem 13.11 and Proposition I3.2[ we are now able to prove our main 
result. 

Proof of Theorem II. 2t First of all, let c G X\ be such that Hcollxf < 1 1 c o 1 1 i? 2 (f2) 

for some fixed constant C > and let 

R = max (C"||co||tf2( n ), \\S(c )L' 1 (cq)J : '(0) \\ Xt ) , 
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where S is defined below. We will construct a solution in Br(0) C Xt for sufficiently 
small < T < 1 with the aid of the contraction mapping principle. If 

ci,c 2 G B R (0) C BUC([0,T};H 2 (Q))nC^{0,T; H\Q)) =: 

< T < 1 with cx\t=o = c-2\t=Qi we conclude 

\\L(d)u - L(c 2 )M|| yT 

< C( ||§( Ci ,©v)-§(c 2 ,Dv)|| l2( o,t;^)+ S(ci,Dv)-S(c2,Dv) 

+ II (7(ci) - 7(c 2 ))vUy,i (ST) + ||(p( Cl )- 2 - p(c 2 )- 2 )Vc , || l2(0iT; h 2 )) 

< C (||ci - C 2 \\ L <x>( Qt T;Bl 1 ) II v II Xi + ll/5(ci)" 2 - p(c 2 ) _2 || L oc( 0iT;/f 2 ) ||Vc' 1^2(0^.^1^) 
+ l|p( c l)~ 2 - P(c2) -2 |Uoo (0iT ._b1 i) ||c / || L 2 (Q)T . H 3 ) ^ 

< C(R) (t^\\ Ci - c 2 \\x2 \\u\\ Xt +T*\\c x - c 2\\Buc([o,n,H^)\\c^\ L4{[o ^ B i^ 

< C{R)T± (J|ci - c 2 || c i ([0)r];ifl) + || c l — C2\\bUC([0,T];H 2 )J \\ u \\x t - 
where u = (v, c') G Xt with c'\ t =o = 0. Here we have used that 



WfgWmm < c (J|/IU 3 x lW IMI^) + \\f\\m { n)\\g\\ Bi(n) 

and have used (13.171) . Hence there is some < Tq < 1 such that 



1 

4(7(1) 



||L(c)u - L(co)u||y T < jttttt II^IUt for all < T < T , ||c|| X 2 < i? 



||L( Cl )u-L(c 2 H| rT < -^—\\u\\ Xt foraU0<r<T ,|| Cj || X 2 <i2,j = l,2, 

since c G X 2 ,. This implies that L(c) : — > is invertible and ||L(c) _1 ||£(y Ti x T ) < 
§C(1) < 2(7(1) as well as 

\\L{ Cl )- x - L{c 2 )- x \\ c{Yt , Xt) 
< 4C(l) 2 ||L( Cl )-L( C2 )|| £(XT , yT) <C(i?)Tl|| Cl -c 2 || x| . 

Moreover, we can choose T so small that llcll v-2 < R and cL=o = Co implies lie — 

T o 

c o|lc°(<2^) < eo since X|> o (7s ([0, To]; (7°(0)), where £o is as in Assumption 11.11 
Then \c(x, t)\ < 1 + e for all < t < T, x G H. 

Hence we can write f l3.7p - fl3.10p as a fixed point equation 

u = 5'(c)L _1 (c) J-"(tt) =: g(u), 
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3 SHORT TIME EXISTENCE OF STRONG SOLUTIONS 



where S(c) : Xt — >■ Xt is defined by 
In order to estimate S'(c), we use 

ll(p( c l) _1 - P(c2) _1 )c'||L2(0,T;H3 (n)3 < C^T^d - C 2 || X 2 || c' \\ X 2 

if ||cj|| X 2 < R, which can be shown using that Xj, C*([0, T]; i?2^ 2 (fi)) and ci — 
c 2 \t=o = c'\ t =o = 0. Moreover, we have 

pmci)- 1 - P{c 2 )- l y)\\ L , {QT) 

< Mci)- 1 - p{c 2 y l w\\ L * iQT) + ii {d t {p{ Cl )- 1 - p( C2 )-y || i2(QT) 

< C(i?) ^||ci - C2|| BC/c .([ 0ir];C .(n))||c / ||x2 + ||9t(ci - c 2)||L 2 (Q T )||c / || jB[ / C .([ 0i T];C(n))) 

< C(R)T±\\ Cl - c 2 || X 2 ||c'|| X 2, 

provided that ||cj||x2 < -R and c\ — C2 1 f=o = c'\ t =o = 0. Here we have used that 

IM Wnc(Ti)) < ^lMlla ([ o,T ]; <(,)) ^ 

for all d G X|. with d\ t=0 = 0. Altogether this implies 

\\S( Cl ) - 5(c 2 )|| £(x o iXt) < C(R)T k q Cl - c 2 \\ X 2 

provided that ||cj|| X 2 < R and C\ — c 2 \t=o = 0, where X^ = {(,c') G X T : c'| t=0 = 0}. 
Therefore we get 

< ll^llr^IlL- 1 ^)-^- 1 ^)!!^^) 

+ \\S{c 1 ) - S(c2)\\c{X^X T )\\L~\c 2 )\\c{Y T ,X T ) 

< C{R)T\\\ Cl -c 2 \\ x p 

where Y$ = {(f , g, a, v , c^) G Y T : Cq = 0}. - Note that L{c)~ l {Y^) = x t- ~ Because 
of Proposition 13.21 we have 

WS^L- 1 ^)^) - S(c 2 )L-\c 2 )F(u 2 )\\ Xt 

< WSi^L- 1 ^) - S{c 2 )L-\c 2 )\\ c{yIXt) \\7{u 1 )\\ Yt + C\\7{ui) - Hu2)\\x T 

< C{R)T*\\ui-u 2 \\ Xt + C{T,R)\\u x -u 2 \\ Xt < Iwm-uzWxr 
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for all sufficiently small < T < T and all Uj = (vj, Cj) G X? with ||ttj||x T < R and 
Cj | i= o = Co. Moreover, 

< \\S(c)L-\c)F(u) - S(co)L-\co)F(p)\\ XT + II^L-^co)^)!!^ 

< ^\\u\\x T + -<R 

for all ||tt||x T < R- Hence by the contraction mapping principle there is a unique 
solution u = (v,c) of flX7j) -f lXTU|) with ||m|| Xt < R. 

Thus we have proved that for every v G H^(Q),cq G H 2 (Q) there are some 
T, R > such that the system fl3.7p - fl3.10p has a unique solution (v, c) G X T with 
||(v,c)||x T < -R- In order to show that there is only one solution (v, c) G X T 
of fl3~7l)-f l3~T0|l . let (v', d) G X T be a second solution of fl3~71)-f l3TT0jl and let R' = 
max(i?, || (v', c')||x T ). Then by the arguments above there is some T" G (0, T] such 
that fl3.7p - fl3.10p has a unique solution (v",c") G Xt> (on the time interval (0,T')) 
with ||(v",c")IU T , < Hence (v,c)| (0 ,t') = (V, J)\ {0 ,t') = (v",c"). Repeating this 
argument finitely many times (with a shift in time), we conclude that (v, c) = (v', c') 
on the full time interval (0, T). ■ 



4 Linearized System — Proof of Theorem 13.1 



In this section we will show unique solvability of the linear system 

d t v - div S(co, Dv) + -|-V div(p" 4 Vc') 

pa 



(n-§(c ,©(P (T v))) r + 7(co)(P CT v) r 
n • v\ 9 n = d n c 



dQ 

an 



fl 


in Qt, 


(4.1) 


h 


in Q T , 


(4.2) 


a 


on S^, 


(4.3) 





on 5t, 


(4.4) 


( v o, c' ) 


in fl, 


(4.5) 



[v,c)\ t=0 

where (fi, f2, a, v , c' Q ) G Yp, (v, c') G X^, and Xt, Yt are as in Section [31 

First of all, we can reduce to the case a = by subtracting from v some w G Xj* 

such that (n ■ S(c , Dw)) T + 7(c )w r = a T , w| 4=0 = 0, divw = 0, wU^ = 0. The 

an 

existence of such a w G X T , depending continuously on a G H^{StY with a n = 0, 
follows directly from Lemma 12.31 For the following let 



T 7 u = (n • S(c , Du)) T + 7(c )u 7 



r)0 



for all u G H 2 {Vt) d . 



Now we will reformulate the system above in an appropriate way assuming that 
(v, c') G Xt- Since (14.21) depends only on div v, we will use the Helmholtz decompo- 
sition to decompose v. More precisely, using L 2 {Vt) d = Ll(Q) © G*2(^) and applying 
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4 LINEARIZED SYSTEM - PROOF OF THEOREM ?? 



P a and (J — P a ) to (14. ip . we obtain that ( 14. ip is equivalent to 

d t w - P a div§(c , Dw) - P a div§(c , V 2 G(divv)) = P a f u (4.6) 

a ^ VG'(divv)-(/-P (7 )div§(co,Dv) + ^Vdiv(p 4 Vc , ) = (/ - P a )f x , (4.7) 

ap 

where w = P CT v and v = w + VG(divv) and G is defined by (I3.4l) - (l3.5p . In the 
following let g = div v. 

In order to determine the principal part of ( 14.61) . we use that 

P CT div§(c , V 2 G(divv)) = P a div (2z>(c )V 2 G(divv)) + P a V(fj(c ) div v) 
= P CT Vdiv((2z>(c ))VG'(divv)) - P a div(2W(c ) <8> VG(divv)) 
= -P a div(2W(c ) ® VG(div v)) = B x g. (4.8) 

Moreover, testing g^J) with V<p, where tp G C£°(0, T\ H l {()) (9)), one sees that (14. 7p is 
equivalent to 



idtg,(p)v, 0T . H -i ) - div§(c ,Dv), 



+4(Vdiv(p 4 Vc),V^) QT = {{I-P a )f x ,V V ) QT 

ap 

for all (p G C^°(0,T; H} Q Jfl)), where we have used again the orthogonal decomposi- 
tion L 2 (n) d = Ll(Q) © G 2 (n) and the fact that g = div„,v G H x (0, T; H^(Q)) if 
v G . Moreover, we note that 

div§(c , V 2 G(divv)) = div(2z>(c )V 2 G(divv)) + V(fj(c )g) 
= 2v{c )Vg + V{fj{co)g) + 2( W(c )) ■ V 2 G(div v) 
= V((2z>(c ) + fj(co))g) + 2(W(c )) • V 2 G(div v) - {2VV{c G ))g. 

Hence 



divS(c Q ,BVG(g)), Vy?J = -(A N (a(c )g), <p) H -i^ + {B 2 g,w) H -^ H y 

where a(co) = 2z>(co) + ?y(co) and B 2 is defined by the equation. 

Therefore we can reformulate (I4.2p - (I4.7I) with a = more abstractly as 

h 

d t u + Au + Bu = | div n (J - P CT )fi | =: / (4.9) 
u\ t =o = I go I =: u (4.10) 
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where u = (c',g,w) T , g = divv , v = w + VG(g ), and 
Au = 



At 



Ai A, \ f(c',gf 

-P a div S(c 0) Dw) / V w 

-P~ 1 Pq 
^A^(divp 4 V-) -A N (a(c ) 




A2W (-div n (divS(co,Dw)))' BU 

and the domains of A, B are defined as 

V{A) = V(B) = {(c',g,w) T :(c',g)eV(A 1 ),weH 2 (n) d nLl(n):T 1 w = 0}, 
V{A X ) = (H 3 (tynH 2 N (ty)xHt 0) {tt). 

We consider A and B as unbounded operators on 

H = Hi x Ll(VL) where Hx = H\tt) x H^(Q). 

For the following analysis it will be crucial that B\ and B 2 are of lower order 
compared to A. More precisely, we have: 

Lemma 4.1 Let s G (|, 1] and Cq G H 2 (Q). Then there are constants C(cq), 
C'(co, s) > such that 

II -01511 £3(0) < C'(s,c )\\g\\ HS{n) , (4.11) 
ll^ll^(0)<^(co)lbll^ (n) (4.12) 

for all g G HLJQ), where B\, B 2 are as above. 
Proof: By the definition of B\ and f !4.8p . we have 

1101511^(0) = ||P ff divS(c 0) lD)VG(5))||£2(o ) 

< C\\Vis\\ H i {n) \\VG(g)\\ H i+s {n) < C(s,c )\\g\\ H ^ n) 

for every s > \ since c G H 2 (tt), A^ 1 : H s (tt) -> H s+2 (tt) for all s G [0, 1] due to 
(ED, and ||/5|| H i < C s ||/||h.|| 5 || h i if 8 > f . 
Finally, i?2 satisfies 

||025||^(o) < 2||(Vz>(co))-V 2 G(^)-Vz>(cokl|£ 2( o) 



< C"||Vc || L 6 (n) ||c/|| L 3 (n) < C(co)ll5ll^ (n) 



due to (TJU). 



Because of the triangle structure of A, it is sufficient to prove that —A\ and 
Po-(divS(co, •)) generate analytic semigroups in order to have the same for —A. 
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Lemma 4.2 Let Ai and Hi be as above. Then —A\ generates a bounded analytic 
semigroup. Moreover, \\Ai(d, g) T \\H\ + |w(c')| is equivalent to WgWm + H c '||h 3 - 

Proof: Let Hq = HL^tt) x Hz}(tt). Then A\ leaves H$ invariant. We first show 
that —Ai\h generates an analytic semigroup. To this end we use that 

-Ak = (_ A \ where A = ^A Jv div(p - 4 V-), J B = -A JV (a(c )-). 

Here V(A) = H 3 (n) n H^(Q) n H^tt), V(B) = H^Q). Without loss of generality 
let (3 = 1. (Otherwise replace A 1 by j3A\.) Because of [T21 Theorem 1.1], — AA Ho 
generates an analytic semigroup on Hq provided that the following conditions are 
satisfied: 

HI A, B are positive self-adjoint operators on HZ-. (fl) with dense domains T>(A), 
D(B). A has a compact resolvent. 

H2 V(A^) = T>{B) and there are constants < p\ < p 2 < oo such that p\A^ < 
B < p 2 A^. 

We note that, if these conditions are satisfied, then e _Al ' H o* is an exponentially 
decreasing semigroup of contractions on H equipped with the norm of D^Az) x 
Hz} (ft). In particular, A\ is invertible. 

Let us verify the conditions above. First of all, A, B are positive and symmetric 
since 



{Au,v) H -x = -— (div(p A Vu),v) L2{n) = — (p 4 Vu,Vv) L2(Q) = {u,Av) H -i 
(Bu',v') H -i = (a(co)u',v') L 2 {n) = (u',Bv') H -i 

for all u, v G T>(A), u',v' G T>(B), where we have used (I2.6p . Moreover, with the 
aid of the Lemma of Lax-Milgram and standard elliptic regularity theory one easily 
shows that A and B are invertible. Hence A, B are self-adjoint. In order to verify 
H2, we use that there are constants cq, Co > such that 

c (V«,VM) L 2(si) < {Au,u) H -i = — (po 4 Vm, Vw) L2(n) < C (yu,Vu) L 2 {n) 

since pZ 4 is bounded above and below, where 

(Vu,Vm) L 2 = -(A N u,u) L 2 = ((-A N ) 2 u,u) H -i = \\A N u\\ 2 H -i. 

(0) "(0) 

Hence co(— Ajv) 2 < A < Co(— A^) 2 in Hz}(Q,). This implies that there are ci, c 2 > 
Cx((-A N )u,u) H -l = Ci||(-A A r)5u||^-i (n) < \\A±u\\ H -i = (A2u,u) H -i 



(0) "-"to)^ (0) '"(0) 

2 

( o)V s v - ' ' '"(0) 



< c 2 \\(-A N )iu\\ 2 H -i {a) = c 2 {(-A N )u,u) H -i 
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because of [23; Chapter I, Corollary 7.1]. Thus ci(-A N ) < A* < c 2 (-A iV ). [231 
Chapter I, Corollary 7.1] also implies V{A^) = HhJQ) = V(—An). Moreover, we 
have that there are c 3 , C4 > such that 

c 3 (-A N u,u) H -i = c 3 (u,u) L 2 < (a(e )~V«)-L 2 = (Bu,u) H -i < c A (-A N u,u) H -i 

(0) (0) (0) 

for all u G HhJQ). Combining this with the previous estimates, we obtain p\Ah < 

B < p 2 A^ for some < p\ < p 2 < 00. 

Hence we have proved H1-H2 and conclude that —A\\h generates a bounded 
analytic semigroup on Hq. Moreover, Hi = Hq © {(0, m) : m G M} and A\(I — Pq) = 
(I - P )A 1 =0. Therefore there is some 5 G (f , it) such that A + A\ : T>(A{) -> Hi 
is invertible for all A G £,$. Moreover, the resolvent estimate ||(A + ) — 1 1 1 < |f[ 
for all A 6 Ej follows from the corresponding estimate for Ai|# . Therefore — Ai 
generates a bounded analytic semigroup on Hi. 

The equivalence of norms follows from the invertibility of Ai|# an d the bound- 
edness of A 1 : (# 3 (fi) PI H%(n)) x HL(Q) -> H v ■ 



Remark 4.3 As shown above there is some 8 > | such that £5 C p(— Ai). In the 
special case that v = z/0,77 = r] , and po = 1 are constant, we have a(co) = 2z/ + r] 
and <r(— Ai) consists of the eigenvalues 

y ap e 

provided that < k < 1, where //„ are the eigenvalues of A^A, cf. [121 Lemma A.l]. 
Note that ^ < 9 < tt and 9 — > | as k — > 0. Hence 5 > | above can be arbitrarily 
close to I in certain situations. 

Because of the triangular structure of A, we conclude from the latter lemma: 

Proposition 4.4 —A generates a bounded analytic semigroup on H . Moreover, 
\\Au\\ H + \\u\\h is equivalent to ||w||#2 + \\g\\m + H c IIh 3 ; where u = (c', g, w) T . 

Proof: Let A 7 v = A 7 (c )v = -P a div§(c , Dv) = -P a div(^(c )©v) for all v G 
£>(A 7 ) = {ue H 2 (tt) d nLl(tt) :T 7 u = 0}. Then -A 7 (c ) : D(A 7 ) C L*(Q) 

generates a bounded analytic C°-semigroup because of Theorem IA. II below. 
Moreover, A 7 is invertible and HAvH^^ is equivalent to ||v||#2( Q ) because of [U 
Lemma 4]. Hence there is some 8 > | such that (A + A)' 1 exists for all A G £5 and 

(x , m-1 _ f (A + A^- 1 -(A + A 1 )- 1 A 2 (A + A 7 )~ 

{A + A) -{ (\+ Ay)' 1 

Using 

|A|||(A + A 7 ) _1 || £ ( L 2) + ||A 7 (A + A 7 ) _1 ||£( L 2) < C s for all A G £, 
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we conclude that 

WA^X + A,)- 1 ^ < C||(A + ^)- 1 f 2 || 7/2(Q) 

< C'\\A 1 {X + A i y l i 2 \\Li{vi) < C'W^Wl 2 ^) 

for all f 2 G L^(fi), A G Sj. Therefore we easily obtain 

C 

||(A + .A) _1 ||h < — uniformly in A G S 5 . 
|A| 

Hence —^4 generates a bounded analytic C°- semigroup. Finally, the equivalence of 
norms can be easily shown using the resolvent identity above for A = 1 and the 
corresponding statement in Lemma [4.21 ■ 

Corollary 4.5 Let Co € H 2 (Q). Then A + B generates an analytic C° -semigroup. 
Moreover, \\(A + B)u\\h + \\u\\b equivalent to ||w||#2 + \\g\\m + ||c||h 3 ; where 
u= (c,g,w) T EV(A). 

Proof: The corollary follows from the fact that for every e > there is a constant 
C £ > such that 

\\Bu\\ H <e\\Au\\ H + C e \\u\\ H 

for all u G T>(A) and a standard result from semigroup theory, cf. e.g. [301 Chapter 
3, Theorem 2.1]. The latter estimate follows from 

\\Bu\\ H < C{co)\\g\\ H%{Q) < C(co)|b||l 2(n) |b||| 1(n) < C(c ) IMl|-i (n) l|sll|i ( Q) 

< C(cb)||^u|||||u||i < + C e (ca)\\u\\ H 

i-t— i 3 

for every e > 0, where u = (c,g,w) and we have used Lemma 14. 1[ H*(Cl) = 
(L 2 (p.),H 1 {q.))s^ 2 , and (EZJ). ■ 

Lemma 4.6 Let A,T>(A), H be defined as above. Then 

(v(A),H) h2 = x Lf 0) (n) x (H\n) d n L*(n)). 

Proof: We only need to show that 

(#(0)(^)' ^(0)(^))^2 = ^(0)(^)> 

(H 3 (n)nH 2 N (n),H\n)) h2 = h 2 n (Q), 

(V(A,),Ll(Q)) h2 = H\Q) d nLl(Q), 

where V(Ay) = {u G # 2 (ft) d PI L 2 (ft) : T 7 u = 0}. The first equality follows 
from i?f 0) (fi) = PoH\tt),L 2 0) (tt) = P L 2 (tt), H^{Q) = PoH- 1 ^) and [35, Sec- 
tion 1.2.4, Theorem]. The second identity is proved using that H S (Q) R if^(O) = 
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(1 - Ajv)- 1 // 1 ^), H^Sl) = (1 - Ajv^H- 1 ^), and H^(Q) = (1 - A^L^O). 
The third identity follows from Lemma [A .21 below. ■ 

Proof of Theorem I3.lt As seen at the beginning of Section H] we can assume 
without loss of generality that a = 0. Let (fi, f 2 , 0, v , c' ) G Y T and / G L 2 (0, T; H) 
be defined as in (14. 9 j) and extend / by zero for t > T. Moreover, let go = divvo, 
w = -P CT v , and let u be as in (I4.10p . 

Applying Theorem 12.11 there is a unique solution u of 

du 

—(t) + Au(t)+Bu(t) = f(t), t>0, 
at 

u(0) = u 

and 

\\{d t u, {A + B)u)\\ L 2 {0tOO . H) < C (||/||i2(o,ooifl) + \\uq\\{v{a),h) 1 ^ 

Using that ||w|| L 2 (0jTo;H) < C(T ) (\\d t u\\ L ^ To . H) + \\u \\ H ) for any fixed < T < oo, 
we obtain that u restricted to (0, T) satisfies 



\\(u,d t U,(A + B)u)\\ L 2 {0tT . H) < C(T ) (\\f\\mO,T;H) + \\U \\{V{A),H) 




uniformly in < T < T and (/, w ). Hence u = (c',g,w) solves (I4.2H - ( I477T) with 
G(divv) = A^g. Therefore (v,c') with v = w + VA^g solves (14. II) -( 1431) . which 
implies 




The estimate of (v, c') G X T , the continuity of L(c ) 1 : Y T — > X T follows from the 
estimate above and the equivalence of norms stated in Corollary 14.51 ■ 



A Stokes Operator with Navier Boundary Condi- 
tions 

In this appendix we summarize some results for the Stokes operator with variable 
viscosity in the case of Navier boundary conditions. More detailed information can 
be found in [TJ Chapter 5]. 
We consider 

A 7 (c): X>(A 7 (c)) C Ll(tt) ->• L 2 a (n): v h-> A 7 (c)v := -P 2 div(2v(c)Dv), 
where c G Wg(fl), q > d, and 

V(Ay(c)) = {u G H 2 (n) d n Ll(Sl) : (2n • p(c)Bu) T + 7 (c)u T = 0} . 
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A STOKES OPERATOR WITH NAVIER BOUNDARY CONDITIONS 



Here we assume that c G W}(fl) for some q > d, 7 G C 1 (M) with < 7(5) < 00 for 
all set, and Q C R d , d = 2, 3, is a bounded domains with C 3 -boundary. - We note 
that H 2 (Q) <^-> Wg(f2) for some g > d if d = 2,3. - If inf sg ]R7(s) = we assume 
additionally that Q does not have an axis of symmetry, i.e., R = {0}, where 

R= { v G H*(n) : v(x) = a + bAx,a,b G M 3 } if d = 3, (A.l) 

fl= |v G H*(Q) : v(x) = a + 6 (^^J ,aG M 2 ,& G eJ if d = 2. (A.2) 

In this case we have the Korn inequality 

||v|| H i < C||Bv|| L 2 for all v G #*(fi), (A.3) 

cf. Necas [221 Theorem 3.5] for the case d = 3. If d = 2, the inequality follows 
from the three dimensional estimate by extending v G H^(Q) 2 to v(xi, x%, X3) = 
(v 1 (a; 1 ,a; 2 ),v 2 (xi,x 2 ),0) T G i#(fi x (-1,1)). 
Because of 

- (div(2i/(c)Dv),w) L2(n) = (2z/(c)Dv,©w) i2(n) + ( 7 (c)v,w) L 2 (an) (A.4) 

= -(v,div(2z/(c)Dw) L 2 (n) 

for all v, w G T>(A 7 (c)), A 7 (c) is a symmetric operator. Moreover, if inf sgK 7(5) = 0, 
then u(s) > Vq > and (IA.3[) implies 

-(div(2i/(c)Dv),v) £3( n ) > c ||v||| 1(n) for all v G P(A 7 ) 
for some c > 0. If 70 := inf seK 7(5) > 0, then one obtains 

-(div(2i/(c)lD)v),v) Z 2 (n) > C(7o) (||Bv|||a ( n) + II v IIl2 ( ^)) > C"(7o)||v|| H i (n ) 

because of ||w||#i(n) < C||Dw|| i 2(Q) for any w G Hq^QY. Furthermore, we have: 

THEOREM A.l Let c G W^(Q), q > d. Then A 7 (c) is a positive self-adjoint 
operator on L^(fi). 

Proof: First of all, A 7 (c): V(A 1 ) -»■ Z£(fi) is invertible because of the following 
arguments: By the Lemma of Lax-Milgram for every f G ^v(^) there is a unique 
v G H^iyt) PI -^v(fi) such that 

(2i/(c)Dv, Dw)i2 (sl) + (7(c) v, w) L 2 (an) = (f , w) L 2( Q ) for all w G T?(A 7 ). 

Since 7(c)v| an G Hh(dti) d , there is some u G # 2 (fi) d n #o( fi ) d n L l(^) such that 
(n • 2u(c)3u) T = — 7(c)v|an due to Lemma [2731 By GauB' theorem we obtain that 
v = v + u solves 

(2u(c)Bv, Dw) L 2 (fi) = (f - div(2z/Du) , w) i2(n) for all w G £>(Ay) . 
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Because of [1, Theorem 5.2.3], we conclude that v G H 2 (Q) d , which yields v G 
H 2 (Q) d . Since Ay(c) is symmetric due to ( 1A.4I) . A 7 (c) is self-adjoint. ■ 

Finally, we need: 

Lemma A. 2 Let 7: R — > [0, 00) 6e continuously differentiable, c G W^fi) /or some 
q > d, and assume that Q possesses no axis of symmetry inf sG R 7(5) = 0. Then 

(Ll(Q),V(A,)) h2 = H\Q)nLl(Q). 

Proof: We use that A^(c) is an invertible, self-adjoint and positive operator on 
L 2 (Q). Hence we can use Theorem 12. II to conclude that for every uo G H 1 (Q)r\L 2 r (Q) 
there is some u G -£^(0, 00; L 2 (fi)) D L 2 (0, 00; T'(v4 7 ) such that u| i=0 = u . Thus 

(Ll{n),v(Aj)^DH\n) d nLl(n) 

by the trace method of real interpolation. But the converse inclusion holds since 
for every u G -£^(0, 00; L 2 D L 2 (0, 00; P(yl 7 )) we obviously have u = u| i=0 G 
H\tt) d n Ll(Q) because of (L 2 (fi), ^(O)) 1 2 = H\n) and (Q. ■ 
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